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Resumen 


MERCADO, J.R., ÍÑIGUEZ, M., GUIDO, P, RAMÍREZ- 
LUNA, J. € GONZÁLEZ-CASILLAS, A. El coeficiente de 
descarga y la densidad beta. Tecnología y Ciencias del Agua. 
Vol. V, núm. 2, marzo-abril de 2014, pp. 161-175. 


Se estudia el coeficiente de descarga y la distribución de 
intensidades de la turbulencia. Con el teorema de Torricelli 
y la teoría de probabilidades se formulan el caudal y el 
coeficiente de descarga, siguiendo una densidad beta 
unimodal, renormalizada, con dos parámetros de forma. 
Se había construido un modelo multifractal para la cascada 
de la energía cinética en la turbulencia, partiendo de los 
métodos de Pearson y de Kolmogorov. Para la intensidad de 
la turbulencia, con el primero se creó una distribución beta; 
para el segundo, una ley en potencia. Se completa el modelo 
multifractal, reconociendo la función de estructura como la 
función Kummer. Se busca la compatibilidad entre los dos 
modelos y la identificación de sus parámetros. Se encuentra 
que los dos parámetros de forma determinan la resolución 
del modelo de cascada. Se determina la dimensión local y el 
espectro de dimensiones para los estados que producen el 
teorema de Torricelli. Redefiniendo la función de estructura, 
la resolución queda determinada por el tirante para el 
cambio de régimen. Análogamente, pueden identificarse 
diversos prototipos, a los que hemos denominado: cuatro 
experimentales, tres canales, Kolmogorov, Kármán, 
Taylor, Verhulst (logística), Cauchy-Manning y Euclides 
(áurea). Se concluye que el coeficiente de descarga es una 
beta renormalizada; la distribución de intensidades de 
la turbulencia es una beta; el prototipo Torricelli resulta 
representativo para los cuatro experimentales y el de 
Euclides, quedando lejos de la distribución Gaussiana, 
que está contenida en el de Kármán; en tanto, el de Taylor 
produce la delta de Dirac. 


Palabras clave: ecuaciones de descarga, autoafinidad, 
modelos de turbulencia, funciones de distribución, función 
de Kummer, multifractales, energía cinética. 


Introducción 


Existen muchos trabajos sobre la determinación 
del coeficiente de descarga para flujos a super- 
ficie libre y en diferentes documentos se pro- 
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Discharge coefficient and turbulence intensity distribution are 
studied. With Torricelli's theorem and the approach of probability 
theory, flowdischargeand dischargecoefficientequation arederivate, 
following an unimodal Beta density function, renormalized, with 
two shape parameters. A multifractals model for the kinetic energy 
cascade in the turbulence was build, starting from the methods of 
Pearson and Kolmogorov. For turbulence intensity, with the first 
method, a Beta distribution was created; with the second, a power 
law. The multifractals model is completed, recognizing the structure 
function as a Kummer function. The compatibility between the two 
models are searched and so the identification of its parameters. It is 
found that the two shape parameters determine the cascade model 
resolution. Local dimension and dimension spectra are determine 
for the two states that produces Torricelli theorem. Redefining 
the structure function, resolution is defined by the water depth 
for the regime change. Analogously, different prototypes could be 
define, which we have call: the four experimentals, three channels, 
Kolmogorov, Kármán, Taylor, Verhulst (logistic), Cauchy- 
Manning, and Euclides (golden proportion). We conclude that the 
discharge coefficient is a renormalized Beta; turbulence intensities 
distribution is a Beta; Torricelli prototype results representative 
for the four experimentals and the Euclides, far away from the 
Gaussian distribution that is contained in von Karman model, 
meanwhile the Taylor 's model yield the Dirac function. 


Keywords: Discharge equations, self-similarity, turbulence 
models, density function, Kummer function, multifractals, kinetic 
energy. 


ponen ecuaciones para flujo libre y ahogado. 
Existe además la deducción clásica de ecuaciones 
de descarga para vertedores de cresta ancha en 
canal rectangular (Henderson, 1966). En Sotelo 
(1999) se citan diferentes ecuaciones de descarga 
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para vertedores por métodos empíricos. En 
dichas ecuaciones se observa la dificultad al 
pasar de un régimen a otro, con coeficientes 
empíricos que conserven su validez cuando 
ocurre dicho cambio de régimen. En general, 
para las estructuras de control y extracción no 
existe una ecuación de descarga que trabaje 
de forma adecuada en todos los regímenes de 
funcionamiento impuestos por la dinámica del 
flujo en los canales de riego. Una ecuación de 
descarga con invarianza para los diferentes 
regímenes es útil para su acoplamiento con 
la dinámica del flujo en el canal. De acuerdo 
con Baume (1992), la derivación de ecuaciones 
de descarga invariantes, bajo los diferentes 
regímenes de funcionamiento, ha sido 
raramente tratada. 

La derivación de ecuaciones de descarga se 
ha abordado con diferentes enfoques empíricos 
y visión determinista, razón por la cual esas 
ecuaciones contienen coeficientes de descarga 
empíricos que dependen de varios parámetros 
(Henderson, 1966; Chow, 1959; Swamee et al., 
1993; Baume, 1992). 

Los regímenes del funcionamiento para el 
flujo en un vertedor se acostumbran denominar 
como flujo ahogado, o a descarga sumergida, 
cuando dicho flujo depende de las condiciones 
hidráulicas aguas abajo y arriba de la estructura; 
en tanto que si sólo depende de las condiciones 
hidráulicas aguas arriba, se denomina flujo 
a descarga libre (Ramírez, 1997). El cambio 
de régimen corresponde al estado de mínima 
energía. Además, el flujo, de acuerdo con la 
estructura, se clasifica de dos maneras: en un 
vertedor, como flujo a superficie libre, y en una 
estructura tipo orificio, como flujo en carga. 

El propósito de este trabajo es estudiar 
el coeficiente de descarga, desde un punto 
de vista experimental y teórico, y abordar 
su invarianza para diferentes vertedores y 
en diferentes regímenes de funcionamiento, 
posibilidad que se abre a través del enfoque 
estadístico. Kármán estudió la función de 
correlación a partir de la ecuación de Navier- 
Stokes; aproximando hasta la doble correlación 
y suprimiendo desde la triple correlación 


(correlación entre tres componentes) en 
adelante y bajo la hipótesis de la turbulencia 
isotrópica llega a una ecuación de Kummer, 
cuya transformada de Laplace inversa produce 
una distribución beta (De Kármán y Howarth, 
1938). Nosotros invocamos las simetrías de 
la ecuación de Navier-Stokes para encontrar 
algunas propiedades de la función de 
correlación y completamos su descripción por 
medio de un proceso de Cantor generalizado, 
(Mercado, 2008). Ahora, con el coeficiente de 
descarga descrito a través de una distribución 
beta, y una transformada de Laplace directa, 
comprobamos la función de correlación, la cual, 
por el trabajo de Kármán, refleja la ecuación de 
Navier-Stokes como una aproximación hasta la 
doble correlación. Además, las simetrías de la 
ecuación de Navier-Stokes son independientes 
de la viscosidad, por lo que la función de 
correlación vale también para la ecuación de 
Navier-Stokes fraccional. 

El esquema del trabajo se presenta en cinco 
apartados. En el primero se define el concepto 
de coeficiente de descarga y se desarrolla el 
método esbozado en esta introducción, la 
cual se sustenta en la ecuación de Navier- 
Stokes, pasa por la ecuación de Bernoulli y 
llega al teorema de Torricelli. En el segundo 
se aplica el teorema anterior y se elabora una 
forma covariante para el caudal adimensional 
generalizado; se reconoce la densidad beta y se 
muestra el cambio de la energía específica de la 
sección, en términos del tirante, hasta su forma 
beta. En el tercero se estudia la relación del 
coeficiente de descarga con la función Kummer, 
conocida como hipergeométrica confluente, y 
su descripción como función Fox; con el modelo 
(Mercado, 2008) se establece la relación entre los 
parámetros de forma de la beta y la resolución, 
para culminar la construcción del multifractal 
Torricelli, identificando la función de estruc- 
tura como función Kummer. En la cuarta 
subsección se estudia la función de correlación 
a partir de las simetrías de la ecuación de 
Navier-Stokes hasta proponer la variable 
dual al tirante y la función de correlación con 
un exponente que contiene la información de 
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los números característicos, y los números de 
Reynolds y Froude. En la última subsección se 
describen los prototipos y los números que los 
identifican, y se construyen tablas que recogen 
esa información; se estudian los números de 
Euler para diferentes vertedores y se identifican 
los valores de sus parámetros, después se 
enuncian las principales conclusiones. 


Métodos 


El coeficiente de descarga puede concebirse 
como un número de Euler que se expresa en 
función de las longitudes relativas carac- 
terísticas y de los números de Froude y 
Reynolds; es un objetivo primordial el escla- 
recer esa dependencia. Los números de Euler 
y de Froude se definen, respectivamente, por: 


E” = y, Ei== (1) 


donde p, v, Ap, 2 y l son, respectivamente, 
densidad, velocidad, variación de presión, gra- 
vedad y una longitud característica. 

En la referencia de Mercado et al. (2012) se 
expone una visión fraccional de la ecuación de 
Navier-Stokes y se considera su versión clásica 
para el presente trabajo, cuya forma vectorial 
es (2): 


0 1 Pp 
—v=vxrotv-V-(-vVv)-V viE+ 2 
rd vVxrotv (-vvv) (5 v e 3 (2) 


donde v, p, p, v, t y q son velocidad, presión, 
densidad, viscosidad cinemática, tiempo y 
potencial escalar de la fuerza externa, como 
la gravedad. Integrando a lo largo de una 
línea de flujo entre dos cualesquiera de sus 
puntos, y bajo la hipótesis de flujo estacionario 
e inviscoso, surge la ecuación de D. Bernoulli 
(1738). Se trata de una ecuación de la energía 
y expresa la constancia de la misma. La 


p 


o (1 mE 
expresión 6 Vivi + > se identifica como 
Po Jia2 


carga total y se establece que la carga total es 
la misma a lo largo de cada línea de corriente 
(Rouse, 1946). La aplicación de la ecuación de 
D. Bernoulli para diferentes estados a lo largo 
de una línea de corriente puede expresarse 
como variaciones entre dos de estos estados, 
cualesquiera de ellos. 

Desde un punto de vista histórico, la ecua- 
ción de D. Bernoulli es una generalización del 
teorema de Torricelli. Este teorema se obtiene 
al considerar que la superficie libre superior 
de un recipiente de agua está en el estado 
que se describe por la presión atmosférica, 
siendo constante y revalorada como p = 0; la 
carga piezométrica por h; y la velocidad v = 0; 
en la otra superficie, la inferior, el estado 
está dado por la presión atmosférica p = 0, la 
carga por h = 0, y la velocidad de salida por 
v. Por tanto, la variación de la energía es 


AE = (0-70) +(0 + gh)-(0+0))=0, entonces 


- e + gh =0, luego v=,/2gh.. Este resultado en 


forma equivalente corresponde a la velocidad 
adquirida por una partícula fluida en el 
proceso de caída libre, que partiendo del reposo 
descienda la altura h bajo la acción única del 
campo gravitacional. 


Coeficiente de descarga 


Para las condiciones del teorema de Torricelli, 
con la caída del tirante H-h, y las diferencias de 
energía a lo largo de las líneas de corriente 
AE =AK + AU, con C como coeficiente de Corio- 


2 
lis, se tiene Ze a o + (h - H) = AE. Para pérdi- 
8 2 
da de carga nula AE=0, entonces H==C hz +h, 


8 


de tal manera que para un ancho unitario, A=h, 


el gasto es Q = AvoQ=h<h EE (m1) la cual 


puede formularse como: 


1/2 
ema 
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Dado que el coeficiente de Coriolis es 
mayor que 1, y puede llegar a acercarse al 
valor 2, su efecto puede interpretarse como 
un reforzamiento del número de Froude, 
visto como una intensificación de la inercia o 
un debilitamiento de la gravedad. En forma 
adimensional, el gasto, para un vertedor de 


cresta ancha en canal de sección rectangular, es 


1/2 
O, =/c-2 l- = ) (Ramírez, 1997). 


22H? A 
h h 


2-1 3/2-1 
ue tiene la forma OQ. =|— e : 
o 2-5) (a) 


Para definir el caudal máximo, según el cambio 


del régimen de funcionamiento, se aplica el 
principio de mínima energía en la ecuación 
adimensional y se obtiene la condición de 
O: 
co 1+1/2 3 
moda), con lo cual se puede definir: 


h 
cambio de régimen 7) (la 


Por lo que la ecuación debe formularse en 
términos de un gasto adimensional, adoptando 
una forma que generalice el gasto derivado 

p a-1 h ba 
del teorema de Torricelli Q. = (5 1-2) , 

H H 
Desde el punto de vista de la teoría de 
probabilidades, el factor (1-H/H)" se ve 
como proporcional al flujo (volumen) por 
unidad de tiempo y por unidad de área y, por 
tanto, representa la densidad de corriente de 
probabilidad; en tanto, el factor (h/H)", como 
proporcional al área, representa la función 
de riesgo, por lo que el producto de los dos 
designa la probabilidad por unidad de tiempo, 
y donde h/H es la cota superior del dominio 
de los valores que puede asumir la variable 
aleatoria en el intervalo (- «o, h/H). Aunque es 
igualmente probable definir al primero como el 
flujo y al segundo como el área, porque bastaría 
con hacer un cambio de variable y medir h/H 


de derecha a izquierda, y no al revés, como 
se acostumbra, lo que equivale al cambio 
h/H>1-(/H); en todo caso, uno de ellos 
representa al flujo de energía y el otro, al riesgo. 


En Ramírez (1997) se formula un coeficiente 
de descarga de la siguiente forma: 


Donde Qu 
descarga libre, y a y b son parámetros de forma, 


es el valor del caudal a 


que dependen de las longitudes relativas 
características y de los números de Froude y 
Reynolds. Si L es la longitud de la cresta del 


vertedor, se observa que Q =C,LO, . . El cambio 
de régimen (+), viene dado por lo BS. 
HA), a+b-2 


la cual corresponde a la moda de la densidad 
de la distribución beta. Por tanto, se representa 
al coeficiente de descarga de forma general 
como: 


Se trata de una densidad beta renormalizada 


Dn” o 
E de tal manera que h = h., valga 1, 
0, 0er), 
por lo que queda normalizada por su máximo 
y no por el área debajo de la curva, como se 
acostumbra con la densidad beta. Ésta es 
la función de distribución más probable en 
el sentido del principio de mínima energía 
y queda determinada por la relación beta 
renormalizada (6), definida hasta el cambio de 
régimen. 

En la sección de exposición de los prototipos 
se ilustra su aplicación al vertedor de cresta 
ancha, de cresta aguda, de compuerta, etcétera. 
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Se conoce el comportamiento característico 
de la energía específica de una sección, que 
muestra un subrango dado por una función 
creciente de la energía con el tirante, el flujo 
subcrítico o lento; otro subrango con una 
función decreciente con el tirante es el flujo 
supercrítico o rápido; un punto que los separa y 
a la vez los une corresponde al tirante crítico, y 
al mínimo de la energía con respecto al tirante. 

El tirante puede adimensionalizarse con el 
tirante crítico, definiendo h/ h,. Con el cambio 

h/h, 

1+1/h, 
de variación desde los reales positivos hasta 


de variable h = se contrae el intervalo 


el intervalo unitario 0<h/h,<w%H0<h<1 
y con el cambio de signo de la energía 
específica de la sección HHH=-H. La 
rama creciente de la función representa el 
subrango supercrítico y la rama decreciente el 
subrango subcrítico, separados por el punto 
máximo, representando el estado crítico. Así, 
el cambio de la densidad de probabilidad es 


a-1 b+1 
eno a ZA e AE 
1+1/h, 1+4/h, 


Coeficiente de descarga y la función 
Kummer 


La representación obtenida en la sección ante- 
rior para el coeficiente de descarga permite otra 
interpretación del mismo, que lo vincula con la 
ecuación diferencial de Kummer. El coeficiente 
de descarga es la transformada de Laplace 
inversa de la función Kummer, en tanto que 
es una solución de la ecuación diferencial 
homónima. O recíprocamente, la transformada 
de Laplace del coeficiente de descarga es 
solución de la ecuación diferencial de Kummer. 
Una sección transversal sobre el vertedor se 
cuadricula y sobre cada cuadro se toma un valor 
del tirante representativo, se adimensionaliza, 
se define la función completamente monótona 
f(1,) s e) ye) con el proceso de Cantor 
(a+b), 
generalizado y la discretización h E —enten- 
diendo como s la variable dual al tirante, la cual 
también debe ser adimensional— se obtiene 


¿EEs; aa + b) (Mercado, 2008). Pero ,F, (2; a,a + b) 
es una solución de la ecuación diferencial 
ordinaria zu”+(a + b — zju' —- au = 0 y además 


Ali” 
B(a,b) 


renormalizada y produce (7): 


¡5 (=s;a,a+b) L 


(s), la cual es 


E(c.(1))( E ¡E (esaa+0)0) 


¡ no (1 > Mo, 


Por tanto, la función ¡le s; a/a + b), solu- 
ción de la ecuación diferencial de Kumwmer, 
determina de forma única al coeficiente de 
descarga, una vez que se han fijado los dos 
parámetros de forma. Sin embargo, de acuerdo 
con nuestra hipótesis, existe una relación entre 
los dos parámetros, relación que depende del 
índice de ocupación espacial o bien del índice 
de estabilidad Lévy. 

Por otra parte, esa función tiene también su 
representación como función Fox y, por tanto, 
como transformada inversa de Mellin (Metzler 
y Klafter, 2000): 


¡E (-s;a,0+b) 


_Flard) na (la) 
T(a) Hz2 (0,1), (1-(a+b),1) 
T(a+b) 


(8) 


r(-1) T(a+t) 
T(a+b+t) 1 
función gamma de Euler. 


en donde x(t)= , y T(t) es la 


Multifractal Torricelli 


En Mercado (2008) se presenta un modelo 
multifractal para la cascada de energía de 
la turbulencia partiendo de los métodos de 
Pearson y de Kolmogorov, simbolizando con q 
la llamada base de la resolución, siendo q > 1; 
además, se tienen f(a.), el espectro multifractal; 
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a, la dimensión fractal local; D,, el máximo del 
espectro multifractal y medida del soporte; e, 
la tasa de transferencia de la energía cinética; 
l,, la longitud que corresponde a la dimensión 
de los vórtices iniciales; y n representando la 
enésima etapa del proceso de fractalización. Se 
encuentra que la intensidad de la turbulencia 
14 
módulo de las variaciones de las velocidades 


entendida como la media cuadrática del 


(como sustituto del valor rms, el cual queda 
indeterminado para las distribuciones Levy, 
salvo la Gaussiana), puede describirse por 


24 Ñ 

(an) (09 (apor con 4, = el, 

_f(a) 
E 


T 
tanto se obtiene (9): 


, y D, es la dimensión topológica, por 


12 
L, Ñ (Ja) _ (el, pe (nt (9) 


distribución de 
la intensidad de turbulencia 


Si se asume que la 
sigue una 
distribución beta h”” (1 — kh), se quiere calcular 
el sustituto del valor rms, para lo cual se 


E 1/2 
requiere evaluar la integral S a n?B dh) ; pero 


su cuadrado es la transformada de Mellin 
de orden 20 + 1 de la densidad beta, recor- 
dando que en tanto distribución se define 
intervalo unitario, 
$, 1"Bdh=M(B)(20+1); en términos de 
los símbolos de Pochhammer, el valor de la 
bei (a). _ T(a+20) P(a+b) 

r(a) T(a+b+20) 


(2+D),, 
tal manera que para una etapa bastante 


como nula fuera del 


de 


elevada del proceso de fractalización, ese valor 


ON 


[0 
A) , por lo 
(a+byo a+b E 


que la intensidad de la turbulencia resulta 


se comporta como 


representada por su valor medio elevado a la 
potencia o. Para que exista compatibilidad, la 
base de la resolución queda definida a partir 


de los dos parámetros de forma g*=-% px 
a+b 


y se denotará como k, más adelante. Además, 
el primer momento de la beta es la tangente 
en cero de la función Kummer, en tanto 
función generatriz de momentos; por tanto, 
si se construye un modelo multifractal, 
en donde la función de estructura sea la 
función Kummer, se puede afirmar que la 
transformada de Laplace inversa del coefi- 
ciente de descarga aporta la función de 
estructura. Para completar la descripción 


del modelo multifractal, se obtiene la dimen- 
_ d d 

sión local U===1(5)= 5 lb (=s;a,a+b)) 

a 

ak ¡E s:a+La+b+1) por lo que el espec- 

tro resulta: 


f(a(s) =1 (=s;a,a+b) 


a 
+5 ,¡E (-s;¡a+1,0+b+1 
Pa A ) (10) 


datos de 
inverso de la resolución es 


En particular, para los 


Torricelli, el 


= e = e entonces la potencia debe ser 
712. 77 P 
o(a)= E - £ -0.19048<1. 
75 =1iaj7 21 


Por tanto, del teorema de 'Torricelli y 
la cascada energética puede construirse la 
distribución de las intensidades de la turbu- 
lencia y renormalizarla para obtener el coefi- 
ciente de descarga, para que entonces su 
transformada de Laplace aporte la función 
generatriz de momentos, la cual es la función 
Kumwmer. Pero, además, para cada uno de los 
prototipos puede construirse el correspondiente 
multifractal con base en los parámetros de 
forma a y b, y la función de Kummer, de una 
manera similar al que hemos denominado 
de Torricelli. Una forma alternativa es definir 
la función de estructura como combinación 
lineal de la variable de reiteración s y la 
función de Kummer desplazada en la unidad 
g(s) = ,F,(=s; a,a + b)— 1, así que g(s) satisface las 
tres condiciones g(0) = g'(0o) = 0, g(00) =-1 y los 
dos coeficientes de la combinación lineal deben 
determinarse para cada caso (Liu et al., 2003). 
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Función de correlación 


El presente camino constituye un recíproco 
al anterior, donde se describe la función de 
correlación de la velocidad y puede obtenerse 
la distribución de probabilidades asociada. 
En general, la función de correlación depende 
de la distancia entre los dos puntos y su 
orientación, pero como la ecuación de Navier- 
Stokes es invariante bajo rotaciones, la función 
de correlación depende sólo de la distancia 
entre los dos puntos y del tiempo (Olver, 1993). 

Por otra parte, dentro de las simetrías 
de la ecuación de Navier-Stokes se incluyen 
las relativas a la escala, las que poseen como 
generadores a: G, = x0, + yd, + 20, + tó, y G,= 
td, 48. =00. 19d, = 2p0,, siendo x, y, z, y t las 


coordenadas cartesianas usuales y el tiempo; 
y u, v, w, y p las coordenadas de la velocidad 
usuales y la presión. Para un canal con el eje 
x en la dirección de la corriente principal, las 


líneas del campo vectorial satisfacen e 00 O 


CIS 


En general, como ya se dijo, la función de 


correlación será función de la distancia y del 
tiempo, pero a través de las transformaciones 
de escala dependerá de las longitudes 
características y de los números de Strouhal, 
Reynolds y Froude. Para el número de 
Strouhal resulta la dependencia de la velocidad 
proporcional, la cual, a través del factor de 
fricción, quedará en términos de los números 
de Reynolds y de Froude. De acuerdo con las 
longitudes características hay que tomar en 
cuenta los dos estados hidráulicos posibles: el 
liso y el rugoso. Cuando la rugosidad queda 
sumergida en la capa límite y se tiene el estado 
completamente liso, la escala de longitud de 


turbulencia está dada por la longitud intrínseca 


as , en tanto, en el estado rugoso está dada 


lo 


x* 


por la longitud asignada a la rugosidad k.. 


Los números de Reynolds locales se dan por 


A para el estado liso y Es E para el 


ES 
Lv Ll v 


estado rugoso. También para el estado rugoso, el 
k 


E 


u 
número de Strouhal se da por ut p u_u 


LL  k Un 


E 


lo 


x 


U e 

, así que para los dos 
0 . 
estados ese número invariante se da por la 


y para el liso 


velocidad proporcional, que a través del 
factor de fricción incluirá la dependencia de 
las longitudes características y los números 
de Reynolds y Froude, dependencia que 
finalmente será transferida a la potencia 
de la velocidad. Por tanto, la función de 
correlación se expresa a través de las variables 


adimensionales L y A , en la subcapa inercial; 
0 , 
así de jo: Kármán, 1938; Chen, 2006). 
l 


, 


o. U, 
Con base en los resultados experimentales, 
para el estado rugoso se admite representar 
la velocidad como distribución en potencia 


(1) 
Yo Vo 


El 
: u t : : 
escala se tiene 2 = a o bien repotenciados: 
u t 
o “0 


-1 o 
u lt 
Us to 


siendo m, > Ll el entero más grande, 


ac uetido: Bb) 


con 0 < O < 1. Con los cambios de 


(11) 


tal que el 0, 


m 
o .po 
con 0 < O < 1l, se modifica la escala, 


m 2 om 
INES 
lo Yo Yo 
ld EA 
TT O 


así que 


entonces 


(m,x-1) 
Ela , por tanto en la ecuación (11) 
lo bo —m¿K -(m,x-1) (Hg 
se tiene [£ pe y entonces 
ha 4 
O= E <1 
m_x-1 


El proceso se basa en los criterios de 
Taylor y Kármán. Sea m(o) =m, el mínimo o 
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admitido en la secuencia (o el máximo o”). 


Bajo el criterio mk — 1 = xk de Kármán, luego 


1 1 
K= 7 a 


a did 
K= , por tanto m= 242 y para Kármán, b 
a+b a 


= 1; por tanto, 4 = , con m > 2. Los valores 


propuestos por Kármán y también por Taylor 


son m, = 5 (De Kármán, 1938), sin embargo, 
esta normalización deja fuera importantes 
resultados reportados en Cheng-Lung (1991), 
pero puede escogerse m > 5 para que queden 
cobijados otros exponentes de rango rugoso, 
como el de Manning. Para ello, en la expresión 


de Kolmogorov a= Ss =D; (1 a B) el término 


ED (1-8 debe tender a cero cuando se 
3T 


suprima la turbulencia y se transforme en 


mono-fractal, con B(s = 0) = > = >, por lo que 
2 2 D : ] 
¿Dr (1 5) 2D, (1 2) >0, si D, > 3. Sin 


embargo, el término z es relativo, ya que puede 


ser absorbido por la resolución para sustituirlo 
por otro, así que la expresión de Kolmogorov 
sea a-=+2(1-p) y la condición D, = 3, en 
cuyo caso m — 2 = 5, por tanto m = 7 es la cota 
inferior y el máximo o* admitido es 6 0 0 = z 


que corresponde al resultado experimental 
propuesto por Manning y está reportado en 
Cheng-Lung (1991). Está 
documentado el vínculo entre la cascada 


suficientemente 


energética de Kolmogorov y la distribución 
Gaussiana para las velocidades (Chen, 2006). 


; : 1 
Este vínculo se mantiene, porque si «(a)=3 


entonces a. = 2.0. Recíprocamente, si se toma 


b =1, para mimetizar el factor (1 — h)'", como 
se hace implícitamente en el modelo de 


: a 1 
cascada, y el inverso de la resolución es k = E 


tal como se propuso en el modelo de cascada. 


Por otra parte, si el máximo aceptado para 


el inverso del exponente de la velocidad es 
a 


O máx 


m y bajo el criterio de Kármán mx -—1 =k, 


luego eE pot tanto m=2+2, y si m=7, 
m-1 a 


dead . b 
criterio de Taylor y Kármán, entonces —<5, 
q 


K 


luego 0< E <5. Entonces, de O= <1, 
q 


m_K-1 


ope 2 a] q 
se infiere que Z=m-0" , y como K=—, 
K a+b 


entonces >> m-o"*-1, Por tanto, bajo el 
criterio Kármám-Taylor: eS 6-07, 

Por ejemplo, si se examina el exponente 
ligado a Lacey k=> conduce a 6 - 4 =2; lo 
satisface. Por el contrario, el de Blasius k = 5 se 
tiene 6-7=-1É (0,5) ; no satisface este criterio 
de clasificación. Por su parte, el de Manning 
K= Ñ produce 6-6=0€ (0,5) y lo satisface 
en el límite. En conclusión, desde un punto 
de vista físico podría señalarse a la fracción - 
como la fracción que separa los exponentes 


en los rangos lisos y rugosos. El intervalo 


experimental reportado por Agroskin> <0O< > 
para el extremo 5 se tiene 6-2 =4 € (0,5), y 
4 7 17 

ara — es 6-—=E/[0,5); por tanto, todo el 
dá 7 4 4 ) E 

intervalo de Agroskin queda incluido y puede 
ser agregado al rango rugoso (Agroskin, 1980). 


Finalmente, la propuesta para la función 
u 

—|=f(s) siendo 
k(9) 


LA Y 
de correlación es f E 


ola 00 


2 
5 E E y o(a) 
rio 


normalizada mediante: si a => 2, se tiene 


o(a) = 1; y para Taylor xk > - > 0(a) >, que 


simboliza una delta de Dirac. 
Resultados y discusión 


En el presente apartado se describen los 
prototipos, los números que los identifican y 
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se construyen los cuadros que muestran esa 
información. Además, se estudian los números 
de Euler para diferentes vertedores y se 
identifican los valores de sus parámetros. 


Prototipos 


Se consideran a continuación varios prototipos 
que provienen de muy diversas fuentes. 

Para ellos, si se aportan los parámetros de 
forma, se obtiene el inverso de la resolución 


dada por e. En seguida puede hallarse 
a+ 


el exponente de la velocidad o(a)= 


ba <1. 
7-1 


Después se encuentra el índice Lévy resolvien- 
1) 

0 
“la 

a 


Finalmente, el índice de ocupación espacial 


do para a. la ecuación k(a) = É - ca) 


(Ramírez et al., 


se obtiene de B(a,D;) =1- 
2009). pes 

Al prototipo principal, los autores lo han 
llamado Torricelli, porque se obtiene del teo- 
rema homónimo. Este prototipo surge de 
los estados que dan origen al teorema de 
Torricelli y además se realiza como el ver- 
tedor de cresta ancha. Los exponentes que 
definen la densidad beta se determinan de 


4 
a =2, b=3/2, luego ia. 


= , en seguida 
a+b 7” su 


b 
E 2 h E 
o(a)=> Ea" y a ==. 
Kk- H),, a+b-2 3 


De la referencia Mercado (2008), reformu- 


2 1 r( 2) 4 
lada como «(0)- (530) z =7 se de- 
QL 


termina que aL = siendo además 


1 
a,D.,)=1- 
p(a,D,) PoR 
Para un vertedor de cresta ancha, y de 


1.1843, 
=0.71854, Dy =3. 


acuerdo con el teorema de Torricelli, si la carga 
total aguas arriba es H, para luego reducirse a 
h, la velocidad del flujo se determina como un 
punto material en caída libre, bajo la exclusiva 
acción de la gravedad, así que v = (2g(H —h))'”, 
y para el caudal en un ancho unitario es 


yapa 3/2-1 
Q=(2gH*) (5 1-5) y en forma adi- 


5 2-1 5 3/2-1 
mensional OQ. =|— 1-= ; en parti- 
Op) (oa) +0» 


sobre la cresta del vertedor será 


OO 
Q 1H H con (+) 


cular, 


2-1 2 
2+3/2-2 3 


H 


El inverso de la resolución es 7 por lo 


que coincide con el mismo valor para el de 
Torricelli; por tanto, el vertedor de cresta 
ancha está representado por el caso que se ha 
denominado Torricelli. 

En el primer grupo se abordan los prototi- 
pos que provienen principalmente de resulta- 
dos experimentales. Unos se sitúan dentro de 
los vertedores, en tanto que otros, como cana- 
les. Para los cuatro experimentales se determi- 
nan los parámetros de los vertedores: cresta 
ancha de arista viva, redondeado y de arista 
redondeada, triangular y Creager (cuadro 1), 
según la ecuación (6), y se anexa el valor de 
las variables al cambio de régimen de funcio- 
namiento hidráulico. Los datos experimentales 
son presentados en Ramírez (1997). 


Cuadro 1. Características de vertedores. 


dores on ansia y anula Cresta ancha a ancha y Eras 
arista viva redondeado arista redondeada 
(a-1) 1.9 1.3 1 0.8 0.618 
(b-1) 0.5 0.5 0.5 0.5 0.5 
(M/H, 0.792 0.722 0.667 0.615 0.553 
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En la figura 1 se muestran las gráficas del 
coeficiente de descarga en el eje vertical y de 
acuerdo con la ecuación (6); en tanto, en el 
eje horizontal se representa (1 — h/H), y 
las gráficas tienen validez hasta el cambio de 
régimen. Se reconocen de izquierda a derecha 
de acuerdo con las columnas del cuadro 1; 
así, la primera línea continua corresponde a la 
primera columna o vertedor de cresta ancha y 
arista viva; la última, de cuadritos, al vertedor 
Creager. En correspondencia, en el cuadro 
2 se muestran los valores de los parámetros 
de forma de la función beta unimodal, la 
resolución, el índice Lévy, el de ocupación 
espacial y el exponente de la velocidad, aunque 


se omite la columna 4 del vertedor de cresta 
ancha redondeada. 

Para los tres canales (Ramírez et al., 2009), 
los datos se presentan en el cuadro 3. 

Los resultados experimentales citados en 
Agroskin (1980), que se ubican en el rango 
rugoso, producen lo que se tiene en el cuadro 4. 

El resultado de Lacey (1930), citado por 
Cheng-Lung (1991), se originó de forma 
experimental, y es una de las potencias más 
altas dentro del rango rugoso. Se obtiene de la 
potencia de la velocidad o = 


==, lo que 
7k-1 4 q 


1 
produce K= 3 , entonces a. = 1.45, por lo que 


B(a,3) = 0.7701. 


Sd 1+ z 0920494 4980 00 
1 ¿eos ór Gob 
+ ++ pa D 
ln] |] 
0.8 
1 + 
+ O 
0.67 
0.4+ 
0.24 
0 + + + h 
0.1 0.2 0.3 0.4 0.5 0.6 17 
H 
Figura 1. Coeficientes de descarga. 
Cuadro 2. Parámetros de los cuatro experimentales. 
Parámetros Resolución Índice Lévy Índice de ocupación Exponente 
1 a=2.9,b= : K = 0.65909 a = 1.1303 P(a,3) = 0.70509 o(a) = 0.18239 
2 a=2.3,b== k = 0.60526 a. =1.1618 p(a,3) = 0.71309 o(a) = 0.18699 
3 a=1.8,b= : K =0.54545 a =1.2031 p(a,3) = 0.72294 o(a) = 0.19355 
4 a=1.618, b 7 k = 0.51892 a =1.224 Pla,3) = 0.72767 o(a) = 0.19713 
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Cuadro 3. Parámetros de los tres canales. 


Resolución Índice de Lévy Índice de ocupación Exponente 

Brannif 0.17845 1.934 0.82765 0.71624 

San Nicolás 0.19413 1.856 0.8204 0.54089 

Coria 0.22114 1.743 0.80876 0.40355 

Cuadro 4. Resultados en Agroskin (1980). 
Exponente Resolución Índice de Lévy Índice de ocupación 

I 1/2 1/5 1.8293 0.81778 
D 4/7 4/21 1.8733 0.82206 


Orificio ahogado 


—pu 5 2 
Para el caso del flujo a través de un orificio, E, = Ap e[r[c.z, | (18) 
siendo Cc el coeficiente de contracción del 
chorro, un poco después de su salida y que está 
b 
2h 


1 
E Al 
E, =h"” (1-h)2” y, en consecuencia, se repiten 


dado por el cociente del área del chorro sobre el 


2 
t 
área del orificio de salida (Rouse, 1946): ] A A 


por tanto, con h-(c 


1 
2P v” C EV los valores del caso anterior. 
E, = =| 2 |=[1-[C, (12) ob 
Ap Co B En el resultado anterior si hb 7 
se transforma en un vertedor de cresta 
por tanto, con los valores de la variable aleatoria aguda, y con 2h>h, el caudal es 
b 3/2 3/2 
ligada al tirante y expresados porh=|C —|,se 2 2 
8 y L Pp C z) q= ¿cda E + ") = 5 (Rouse, 
tiene E, =W” (1- h)2 >. 3 28 28 


En consecuencia, el inverso de la resolución 2 v ús 
1.2 1946), desplazando q/=93+ Ea ze L 
es k(a) DE luego produce a. = 1.1262 3 28 
1+2 2 AD 2 
* es cr , 5-22, o y 
con o(a)=7 y Blo,3) = 0.70402. E 0 A 


donde números F, bajos implican grandes de- 
Orificio horizontal ahogado 
flexiones gravitacionales: 


Para la descarga a través de un orificio 
horizontal ahogado se debe sustituir la 


E 3/2 
dd bobo. - E == =(1 h) (14) 
proporción entre fronteras E > > siendo 2h 0 
relativamente grande con respecto a b, o sea la 
itució h ; 52.3 5/7 5 
sustitución del ancho B > 2h, entonces (Rouse, así que k= e E AñO 
1946): 5/2+1 7 7(5/7)-1 28 
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y se ubica en el rango rugoso y cerca del de 
Manning. Además, se obtiene a = 1.1022; y 
pla,3) = 0.69757. 

Se recuerda que en este vertedor se presenta 
un decrecimiento de la presión piezométrica 
desde el dato, ubicado en la base plantilla, hasta 
el borde del vertedor de altura 1w, en donde se 
hace cero, acompañado con un crecimiento de 
la carga de velocidad; este borde representa 
un cambio brusco, un quiebre de la pendiente 
de la presión piezométrica. Luego dentro del 
chorro, un posterior crecimiento del tirante 
de presión hasta un máximo en el interior del 
chorro, para luego presentar un decrecimiento 
hasta la superficie exterior superior del chorro, 
que vuelve a ser cero, con el concomitante 
decrecimiento de la carga de velocidad, seguida 
de un crecimiento hasta la superficie exterior 
del borde superior del chorro (Rouse, 1946). 


Compuerta 


En este prototipo se presenta una diferencia 
con el de la abertura horizontal ahogada de- 
bido a que la presión en el exterior del chorro 
de aquel es la atmosférica, en tanto que para el 
presente se tiene una distribución hidrostática. 
Con la ecuación de continuidad ¿ la de D. 


Pa 


Bernoulli 
(Rouse, 1946), luego: 


se encuentra q= E 


11e a se reproduce el 
11241 3% "P 


de Lacey. Se presenta un comportamiento del 


Entonces K= 


tirante de presión y de la velocidad similar al 
vertedor anterior con el cambio de la altura w 
por la abertura b. 

El segundo grupo de prototipos proviene 
de la teoría (algunos a partir de los parámetros 
de forma), otros del índice Lévy y algunos de 
la potencia. 


Kolmogorov 


Este prototipo surge de la cascada energética 
I, (ecuación (9)), produciendo un modelo en 
potencia, que corresponde a una densidad beta 
B = h*(1 — hY”, con b = 1 y siendo el índice 
Lévy a = 2, que designa a la distribución 

a 1 
a+b 6 


se obtiene a = > Además, p(2,3) = > En tanto, 


Gaussiana, por lo que k(a) = Ñ y dex= 


ll 
7-1 
régimen 


el exponente de la velocidad es o(a) = 


El tirante cambio de 


h = a-1 =1. El modelo de cascada 
er A+b-2 


para 


se ilustra en la figura 2. En el eje vertical se 
representa la intensidad de la turbulencia 1, y 


en el horizontal h = q”. La curva de trazos es 
1 


a! 
dl , ya que el coeficiente de la función beta 


OS 


Sobre esta curva de trazos, en la secuencia 
de pequeños círculos, se aprecia el proceso 
de fractalización que ocurre de derecha 
a izquierda, en donde la sucesión de la 


n 
AN 1 
resolución se toma como h HP (5 de acuerdo 


con el planteamiento original de Kolmogorov, 
aunque también puede tomarse por ejemplo 


EN ) 
como ho 2 . En tanto, la línea continua 


1 


dl 
representa 1 po y sirve de referencia 
B(1/10) 


para el cambio del parámetro de forma a -- 


pl de : 
hasta a=—, y también para el cambio en la 


SS mm 1Y 
sucesión de la resolución desde h > 5) hasta 


hos E 
10 
El de Kármán surge de la función de 


correlación, al establecer el inverso de la 


resolución como k = E así que en la distribución 


b = 51. La potencia de la velocidad o(a) = 1, 
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0 0.2 0.4 0.6 0:8 1 


Figura 2. Modelo de cascada. 


por lo tanto coincide con la de Kolmogorov. 
El cociente de funciones gammas puede 
expresarse como una función Pochhammer 


y ésta se reduce a la identidad 


si —-—=1, lo que equivale a a. =2, y entonces 
2 1 15 1 1 
lla ALL [1] =yk=2:2=- Defor- 
356 Ql (1) a 


ma recíproca, cerca del régimen laminar a. > 2, 


1/1 LT: 1 
por lo que k= =="===, lo cual 
312 3) 32 6 
22 


significa la participación de la viscosidad. 


Taylor 


El modelo de Taylor brota de la función de 
pas E Laa 
correlación, siendo k(a) == el inverso de la 


resolución. De ello se obtiene a. = 2.1563, lo cual 
no es posible, porque a. < 2 y resulta a = oo, 


Manning 


Si el índice Lévy es a = 1, la distribución es la 
de Cauchy, luego el inverso de la resolución 


es k(1) = 1, entonces a | y b= 0; la 
a+b 
potencia de la velocidad es o-- y el índice 


de ocupación espacial es B(1,3) = = De forma 


análoga, de acuerdo con los resultados de 


Manmning, la potencia de la velocidad es O = A 
lo que produce k(a) = 1, entonces se obtiene 
a=s1ly B(1,5)-5 Este exponente de modo 


tradicional se ubica en el rango rugoso. 
Euclides 


Se origina en el arte y puede describirse a partir 
de la distribución, recordando la definición de 
la proporción de oro de Euclides: un segmento 
puede partirse en dos subsegmentos de 
longitud; el mayor es la media proporcional 


entre la totalidad y el menor da > lo que da 
q 


E 1 1 
lugar a la solución a = (7-34) o en forma 
aproximada surge de dividir un segmento 
en cinco partes y destacar tres de ellas contra 
las dos restantes. Entonces, el inverso de la 


2 -0.61803= > Elindice 


1+ 4/5 


Lévy es tal que a = 1.1652. La potencia de la 
velocidad es o(a)=0.1858 = z El índice de 


resolución es k = 


ocupación espacial es f3 = 0.71393. 
Logística 


Puede originarse desde cualquiera de dos as- 
pectos: las funciones de distribución o de 
correlación. Desde la distribución, sia=2=b, 


entonces el inverso de la resolución es k-5 
y el índice Lévy a = 1.2401. La potencia de la 
velocidad es o(u)== y B(a,3) = 0.7312. En la 
figura 3 se muestra la gráfica de los inversos 
de la resolución con respecto al índice de 
estabilidad, donde se aprecian sus ubicaciones 


relativas y sus posiciones en referencia a los dos 
extremos marcados por a =1 y a. =2. En el eje 
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horizontal se representan los valores del índice 
Lévy a. La línea continua describe la gráfica de 


3 
0-39 a 


"la 


el caso de Taylor, que se ubica fuera del intervalo 


. En particular, se destaca 


señalado. Tanto los cuatro experimentales 
como el de Euclides de la proporción dorada 
se agrupan alrededor de la posición del de 
Torricelli y por tanto se asemejan al vertedor 
de cresta ancha; todos ellos, a su vez, se ubican 
más cerca del extremo 1 de la figura 3, por lo 
que se destaca su comportamiento turbulento. 

La potencia de la velocidad se ilustra en la 
gráfica de la figura 4, en donde se señalan los 
dos extremos representados por a = 1, caso de 


Cauchy h 0(0) = a > > y a. =2, caso de Gauss 
> o(a) =1. En particular, para el de Taylor, la 
potencia diverge a infinito, así que está fuera 
del intervalo 0 < O < 1, como ya se señaló. La 
k(a 

MO 
7x(a) -1 

el eje horizontal se representa el índice Lévy. 


línea continua representa o(a) = 


La ecuación (6) permite representar de 
forma general las ecuaciones de descarga para 
diferentes vertedores. Se conjetura que dicha 
ecuación es válida para otros vertedores, lo que 
se ilustra encontrando la forma de los números 
de Euler para el orificio ahogado, el orificio 
horizontal ahogado y la compuerta. Los datos 


Figura 3. El inverso de la resolución k(a). 
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Figura 4. El exponente de la velocidad o(a). 


de los parámetros de los tres canales muestran 
la validez de las relaciones presentadas en el 
primer párrafo de la sección Prototipos. Los 
parámetros correspondientes al modelo teórico 
Torricelli son representativos de los cuatro 
experimentales, los que además se ubican 
claramente en régimen turbulento (figura 3), 
mientras que, por el contrario, el modelo teórico 
de Kolmogorov resulta más representativo 
para los tres canales, que se ubican más cerca 
del régimen laminar. 


Conclusiones 


1. El coeficiente de descarga adopta la forma 
covariante que se deriva del teorema de 
Torricelli o vertedor de cresta ancha. 

2. La distribución de la intensidad de la 
turbulencia es la beta. 

3. La función de correlación es la función 
Kummer. 

4. Los exponentes de la velocidad del rango 
rugoso se extienden desde el de Maming 
hasta el de Kolmogorov o Kármán. 

5. El coeficiente de descarga (6) refleja, a 
través de una aproximación, las ecuaciones 
de Navier-Stokes. 
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